LARGE DEVIATIONS FOR INTERACTING BESSEL-LIKE PROCESSES 
AND APPLICATIONS TO SYSTEMIC RISK 



TOMOYUKI ICHIBA AND MYKHAYLO SHKOLNIKOV 



Abstract. We establish a process level large deviation principle for systems of interacting 
Bessel-like diffusion processes. By establishing weak uniqueness for the limiting non-local 
SDE of McKean-Vlasov type, we conclude that the latter describes the process level hydro- 
dynamic limit of such systems and obtain a propagation of chaos result. This is the first 
instance of results of this type in the context of interacting diffusion processes under explicit 
assumptions on the coefficients, where the diffusion coefficients are allowed to be both non- 
Lipschitz and degenerate. In the second part of the paper, we explain how systems of this 
type naturally arise in the study of stability of the interbank lending system and describe 
some financial implications of our results. 



1. Introduction 

Our main objects of study are systems of interacting diffusion processes on R + = [0, oo) 
given by 

(1.1) dXi(t) = b(X i (t),p n (t))dt + a(X i (t))dW i (t), z = l,2,...,n, 

where p n (t) = - Y^=i^ x i{t) * s the empirical measure of the system at time t and b, a are 
functions with values in R, R + , respectively. We will be interested in situations, where 
cr(0) = and a is non-Lipschitz around 0. The main example we have in mind is o(x) = 
Cyfx (see the applications below), but our analysis applies to suitable perturbations of such 
functions as well. 

To state our main results, we need the following set of notations. We write Mi (R + ) for 
the space of probability measures on R + endowed with the topology of weak convergence 
of measures, and let Mi(R+) be the space of probability measure on R + with finite first 
moments endowed with the 1-Wasserstein distance. The latter is defined by 



(1.2) d^p, v) = inf / \x-y\dp(x,y), 

where the infimum is taken over all probability measures p on (R+) 2 with first marginal given 
by p and the second marginal given by v. Next, we fix a T > (which will remain fixed 
throughout) and let y = C([0, T], Mi(R + )) be the space of continuous functions from [0, T] to 
Mi(R + ) equipped with the topology of uniform convergence. We set X = M^CQO, T], R + )) 
for the set of stochastic processes with continuous paths taking nonnegative values such that 
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all their one-dimensional marginals have finite first moments. For an element 7 G X, we 
write (7r(7))(t) = ^(t) G M 1 (M + ), t G [0,T] for the path of one- dimensional marginals of 7. 
Finally, we define the topology on X as the collection of preimages of open sets in y, which 
are open in the topology of weak convergence on X (with C([0, T], R + ) being endowed with 
the topology of uniform convergence). In other words, a sequence 7 n , n G N in X converges 
to a limit 7 G X iff 7 n — > 7 weakly and vr(7„) — > 71(7) in y. 

Our main results (large deviation principle, law of large numbers and propagation of chaos) 
hold under the following assumption. 

Assumption 1. The function b is Lipschitz on IR + x Mi(R + ) and it holds 6(0, •) > 0; 
whereas the function a is of the form o~(x) = \fxg(x) with a continuous bounded function 
g : M.+ — > M + taking strictly positive values and satisfies 

(1.3) \a(x) - a{y)\ < 0(|x - y\), x,y G R+, 

where : M + — > R + is such that Q ^ 2 da = 00 for all e > 0. Moreover, the sequence 
p n {0), fi6N is deterministic and converges in Ml(R + ) to a limit A G Mi(R as n — > +00. 

The following are the three main results of our paper. 

Theorem 1 (Large deviation principle). Under Assumption [U the sequence p n , n G N 
satisfies a large deviations principle on X with scale n and good rate function 

u(t) 2 dt , 

where E 1 is the set of (u, W) such that W is a standard Brownian motion, u is a progressively 
measurable process with respect to the filtration generated by W satisfying E J Q T u(t) 2 dt < 
00 and there is a weak solution of 

(1.5) dX{t) = (b(X(t), £{X{t))) + u(t)a(X(t))) dt + a{X{t)) dW(t) 

with law 7. Here, C(X{t)) is the law of X{t) and we use the convention that any infimum 
over an empty set is given by +00. 

As we will see below, Theorem [T] implies the following law of large numbers and propagation 
of chaos result. 

Corollary 2 (Law of large numbers). Under Assumption^ the X -valued sequence p n , n G N 
converges in distribution to the law of the unique strong solution of the nonlocal SDE 

(1.6) dX(t) = b(X(t), C(X(t))) dt + a(X(t)) dW(t), 
where £(X(t)) is the law of X(t) and W is a standard Brownian motion. 

Corollary 3 (Propagation of chaos). Suppose that Assumption^ holds and that for some k G 
N the initial conditions (Xi(0), X 2 (0), . . . ,Xfc(0)) converge to a limit (x\,X2, ■ ■ ■ ,£&)• Then, 
the law of the process (X 1; X 2 , . . . , X^) converges to the product measure fi Xl ®p X2 ®- ■ -®p Xk , 
where p Xi , i = 1, 2, . . . , k are the laws of strong solutions to 

(1.7) dXi(t) = b(Xi(t), p°°(t)) dt + a(Xi(t)) dWi(t), 

started from Xi, i = 1, 2, . . . , k, respectively, and p°° is the law of the solution to ( 11. 61) with 
initial condition distributed according to A. 



1.4) J( 7 ) = - inf E 

' 2(u,W)6B 7 
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Interacting diffusions, for which the drift and diffusion coefficients of each process depend 
only the current value of that process and the current value of the empirical measure of 
the whole system, are usually referred to as diffusions interacting through their mean-field. 
The behavior of such processes has been a classical object of study in the statistical physics 
literature (see [7] and the references there) following the seminal work [11]. However, large 
deviation results for this kind of interacting particle systems are rarely found. Notable 
exceptions are the articles [3J, [I] and [2J. 

In [3J, a large deviation principle is established for interacting diffusion processes as in (II. ip 
under the assumptions of non- degeneracy of the diffusion coefficients and the well-posedness 
of the corresponding martingale problem. In [1] , a large deviation principle is proven for sys- 
tems of diffusions interacting both through the drift and through the diffusion coefficients, 
however for a very special type of interaction known as interaction through the ranks; there, 
also the non-degeneracy of the diffusion coefficients is one of the assumptions. Finally, in [2J 
the authors give a large deviation principle for a very general class of diffusions interacting 
through their mean-field, but under non-explicit assumptions such as the strong existence 
and uniqueness of solutions to ( 11. ip . weak uniqueness of solutions to (11 .5p and tightness of 
empirical measures corresponding to tilted versions of ( II. ip (the latter is equivalent to expo- 
nential tightness of the original particle systems). The latter three assumptions are difficult 
to check in general if the drift and diffusion coefficients are non-Lipschitz or unbounded. In 
our setting, the diffusion coefficient can be non-Lipschitz and unbounded. Nonetheless, we 
are able to show that the suitably modified abstract assumptions of [2J hold in our setting 
and obtain Theorem [1] by adapting the proof in [2J correspondingly. Another main contri- 
bution of the present paper is the proof of strong uniqueness for the non-local SDE (11.61) . 
which allows us to deduce Corollaries [2] and |3J from Theorem [TJ 

In [B] interacting diffusion processes have been proposed as models for the interbank 
lending system with the goal of analyzing the stability or instability of the latter. Following 
this program, we study the features of the weak solution of (11.11) viewed as a model of the 
interbank lending system in the second part of the paper. In this context, Xi(t) stands for 
the monetary reserve of bank i at time t. Banks lend and borrow money from each other in 
order to hold the necessary amounts of monetary reserves. Considering appropriate choices 
of the functions b, a, we illustrate the effect of the monetary flows, which are represented by 
the interactions through the drift coefficients. 

As an application of Theorem [1] and Corollary [2J we shall derive the dynamics of the 
limiting distribution of banks' monetary reserves as n — > oo. We are interested in the effects 
of banks' interactions on (1) the long term behavior of the large interbank lending system, (2) 
the volatility of the limiting system and (3) the boundary behavior of the limiting system 
when most banks are almost bankrupt. We expect that such probabilistic studies of the 
interbank lending system will help understand the key features of systemic risk and some of 
the reasons for financial crises. 

The rest of the paper is organized as follows. In section 2.1 we prove weak uniqueness for 
the limiting non-local SDE ( II. 6p (despite the fact that the diffusion coefficients in (II. ip can be 
non-Lipschitz and degenerate) and thereby derive Corollary |2j from Theorem[TJ Subsequently, 
we show how Corollary [3J can be obtained from Corollary [2J In section 2.2 we prove that the 
suitably modified abstract assumptions of [2] apply in our setting and use the correspondingly 
modified arguments from [2J to complete the proof of Theorem [TJ The main challenge to 
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overcome here is the unboundedness of the coefficients in (11.11) . In section 3 we apply these 
results to model the large interbank lending system by choosing a particular form of the 
functions b, a. The solution of the non-local SDE ( II. 6p then becomes a non-local squared 
radial Ornstein-Uhlenbeck process (see Corollary [6]) . In section 3.1 we discuss the moments 
of the latter, recurrence and transience (see Propositions EHH]), whereas in section 3.2 we 
study the effect of the interactions on the volatility of the limiting system. Section 3.3 is 
devoted to the derivation of a quasilinear first-order PDE satisfied by the Laplace transforms 
of one-dimensional distributions of non-local squared radial Ornstein-Uhlenbeck processes. 
Lastly, in section 3.4 we compute the stationary distribution of the latter, which corresponds 
to the long-term behavior of the limiting interbank lending system. 

2. Large deviation principle 

2.1. Proof of Corollary [2J The proof of Corollary [2] relies on the following uniqueness 
statement for the limiting non-local SDE, which is of independent interest. 

Proposition 4. Let b, a be as in Assumption^ Then, the solution of the non-local stochastic 
differential equation 

(2.1) dX{t) = b{X{t), C(X(t))) dt + a(X(t)) dW{t) 

is pathwise unique. In particular, the law of the pair (X, W) is uniquely determined. 

Proof. Let X, X be two solutions of ( 12. ip with respect to the same Brownian motion W and 
with the same initial value. Define the stopping times 

(2.2) t l : = inf{t > : max(X (t) , X (t)) > L}, LeN. 

In view of Corollary IX.3.4 in [T2], the local time of X — X at zero is identically equal to 

zero. Therefore, Tanaka's formula implies that the process 

(2.3) 

/•'ATI, 

\X(- A t l ) - X(. A Tl)\ — I sgn(X( S ) - X{s)) (b(X(s), C{X{s))) - b(X(s), C(X(s)))) ds 

Jo 

is a martingale starting at zero, where sgn(x) = l{ x >o} ~~ l{x<o}- Evaluating at a fixed t > 0, 
taking the expectation and using the fact that the function b is Lipschitz (see Assumption 
[T]), we obtain 

l>tAT L _ 

E[\X(t A r L ) - X(t A t l )\\ < CE / [\X(s) - X(s)\ + d.mxis)), C(X(s)))] ds 

l Jo J 

with a uniform constant C < oo. Moreover, the definition of the 1-Wasserstein distance 

implies 

d 1 {jr(X(s)),jr(X{s))) < E[\X{s) - X(s)\], s > 0. 
Passing to the limit L — > oo, using Fatou's Lemma for the left-hand side and the Monotone 
Convergence Theorem for the right-hand side, we obtain 

E[\X(t)-X(t)\]<2C [ E[\X(s)-X(s)\]ds. 

Jo 

Therefore, Gronwall's Lemma shows that the function 1 1— > E[\X(t) — X(t)\] vanishes identi- 
cally. This yields the desired pathwise uniqueness. Finally, in view of the theory of Yamada 
and Watanabe (see Proposition 5.3.20 in [10]), pathwise uniqueness implies that the law of 
the pair (X, W) is uniquely determined. □ 
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We are now ready to give the proof of Corollary [21 taking Theorem [T] for granted. 

Proof of Corollary [2l The formula for the rate function in Theorem [1] shows that J{^) = 
holds if and only if 7 is the law of a solution to the nonlocal SDE (12. ip . By Proposition [4] 
this law is unique, and Corollary [2] now follows from the goodness of the rate function. □ 

Next, we show that Corollary [3] is a consequence of Corollary [2j 

Proof of Corollary® We fix k and x±, x 2 , . . . , x k as in the statement of Corollary |3j Corollary 
[2] implies that the sequence of random variables p n , n G N is tight on Mi(C([0, T], 
Therefore, by Proposition 2.2 in [T5] the sequence of processes (X l7 X 2 , . . . ,X k ) as n varies 
is tight as well. Moreover, the same argument as in the proof of Corollary [2] shows that 
the solution of the SDE ( I1.7P is pathwise unique. Therefore, it suffices to identify the limit 
points of the sequence of (Xi, X 2 , . . . , Xk) as n varies as the unique solution of the following 
martingale problem: for each / G C^°(M. k ) the process 



k 



(2.4) /(*(•)) -f(x)- ^ J q b(X t (t),p-(t))^-(X(t)) + -a(X l (t)Y^(X(t))dt 

is a martingale, with x = (xi,X2, ■ ■ ■ ,x k ). To this end, it suffices to note that by Corollary 
[2] each process as in (12 .4p is the limit of the processes 

M^(-) := /(Xx(-), • • • , X k (-)) - /(X!(0), . . . , X k (0)) 



J2 J °(Mt))^(Xi(t), x k (t)) dWi(t) 



as n — > 00, and that sequence of random variables M n, f(t), n G N is uniformly integrable for 
any fixed t G [0, T], since the second moments of the latter are uniformly bounded by virtue 
of Ito's isometry, the continuity of a and the fact that / has compact support. □ 

2.2. Proof of Theorem [J3 This subsection is devoted to the proof of Theorem [TJ We first 
recall the framework of [2] as our approach will rely on a version of the main result there. 

The large deviations principle in [2] is established in form of a Laplace principle, with 
the local large deviations upper and lower bounds obtained by considering suitable tilted 
versions of the original system of SDEs. In our case, the appropriate controlled version of 
the system (II .ip is given by 

(2.5) dX t (t) = (b(Xi(t),^(t)) + Ul (t)a(X t (t))) dt + <r(Xi(t)) dW^t), i = 1, 2, . . . , n. 

Hereby, the control processes u\, u 2 , • • • , u n are assumed to be progressively measurable with 
respect to the filtration generated by the Brownian motions Wi, W 2 , ■ ■ ■ , W n and to satisfy 



< 00, 



pi 

(2.6) V E / Ui{tfdt 

i=i lJ o 

and = ^ Xir=i ^xJt) i s ^ ne empirical measure of the controlled system at time t. 
In addition, we introduce a controlled version of the limiting non-local SDE (II. 6ft : 

(2.7) dX(t) = (b(X(t), C(X(t))) + u(t)a(X(t)))dt + a(X(t)) dW(t), 
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where u is a progressively measurable process with respect to the filtration generated by the 
Brownian motion W, which satisfies 



(2.8) 

and £(X(t)) is the law of X(t). 



E 



u(t) 2 dt 



< oo, 



The exponential tightness of the original sequence p n , n £ N of empirical measures cor- 
responds in this setup to the precompactness of the sets {ff 1 : n E N} for all sequences of 
controls satisfying 



(2.9) 



i n r f 1 

sup - V E / 

neN n ^ [Jo 



\Ui(t)\ 2 dt 



< oo. 



In this situation, a slightly modified version of the main result in [2] reads as follows. 

Proposition 5. Assume the following. 

(i) The initial conditions p n {0), n £ N are deterministic and converge in M\(R + ) to a 
limit A. 

(ii) The coefficients b, a are continuous on R + x Mi(R), M+, respectively. 
(Hi) Strong existence and uniqueness holds for (11.11) . 

(iv) Weak uniqueness in X holds for (12. 7p for any given joint distribution of (u, W). 
(v) For any sequence of controls satisfying (12. 9p . the set {p n : n £ N} is precompact in X . 

Then, the sequence p n , n £ N satisfies a large deviations principle on X with scale n and a 
good rate function 

[ u{tfdt , 
Jo 

where £> 7 is the set of all (u, W) such that the corresponding weak solution X of (12. 7p has 
law 7. 

Proof. The proposition is obtained by carrying out the proof of Theorem 3.1 in [2] replacing 
their space X by our space X. □ 

One of the main contributions of our paper is to show that the abstract assumptions 
(iii)-(v) in Proposition are satisfied under the concrete Assumption [T] on the coefficients in 
(11.11) . so that Theorem [T] becomes a consequence of Proposition 

Proof of Theorem^ In view of Proposition [5], it suffices to check that the assumptions (i)-(v) 
there are consequences of Assumption [lj The assertions (i) and (ii) are immediate, so that 
the rest of the proof is devoted to the derivation of assertions (iii), (iv) and (v), which are 
dealt with in steps 1, 2 and 3, respectively. 

Step 1. To prove (iii), we note first that a weak solution of (II .ip exists by virtue of Theorem 
A in pp. Therefore, it suffices to prove the pathwise uniqueness for the solutions of (II. ip . To 
this end, we let (Xi, X 2 , . . . , X n ) and (Xx, X 2 , . . . , X n ) be two solutions of (ll.ip on the same 
probability space and with respect to the same Brownian motions Wx, W 2 , ■ ■ ■ , W n . Then, 



(2.10) 



J(i) = - inf E 
2 («,K0eB 7 
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by Corollary IX.3.4 in [12] . the local times of X 1 — X±, X 2 — X 2 , . . . , X n — X n at zero all 
vanish identically. Therefore, for any fixed i G {1,2, ... ,n} and LgN, the process 



. (n) 



A rf) - X t (- A r£ n) )| - / sgn(X,( S ) - X t (s)) , p n (s)) - b^s), p n (s))) ds 

Jo 

is a martingale starting at zero, where 

rjp = inf{t > : max maxpQ^), X^t)) > L}. 

8=1,2,. ..,11 

Since the function b is assumed to be Lipschitz (see Assumption [1]) and one has 



n 1 n ~\ n 



j=i j=i j=i 

it follows 

n „. n 

(2.11) ^ E[|X,(. A ri n) ) - A rf 5 )|] < C / J]E[|X,( S A r^) - X,( S A )|] ds 
<=i Jo i=i 

with a uniform constant C < oo. Applying Gronwall's Lemma first and then taking the limit 
L — > oo, we conclude Xi(t) — Xi(t), i — 1,2, ... ,n for any fixed t £ [0, T] with probability 1. 
Thus, due to the path continuity of both solutions it must hold Xi(-) = Xi(-), i — 1,2, ... ,n 
with probability 1. 

Step 2. We now turn to the proof of assertion (iv) in Proposition [5j By the theory of Yamada 
and Watanabe (see Proposition 5.3.20 of [10]) it suffices to show that the solution of (12. 7p is 
pathwise unique. To this end, let X 1 , X 2 be two solutions of (j2.7p on the same probability 
space and with respect to the same pair (u, W). Now, introduce the stopping times 

(2.12) a M := inf {t > : J u{s) 2 ds > m}, M e N, 

and for each fixed M G N define a change of measure according to the density 

(2.13) — = exp ( - ^ u(s) dW(s) --J u{sf ds) . 

In view of Novikov's condition, P M is a well-defined probability measure, which is equivalent 
to P. Moreover, by Girsanov's Theorem it holds 

/••Acta/ /"Acta/ 

(2.14) X 1 (-Aa M )=X 1 (0)+ / b(X 1 (t),C P (X\t)))dt+ / a{X x (t)) dW M (t), 

Jo Jo 

/"Acta/ /"-Acta/ 

(2.15) X 2 (-Aa M ) = X 2 (0)+ / 6(X 2 (t), £ P (X 2 (t))) dt + / a(X 2 (t)) dW M (t) 

Jo Jo 

under P M , where C P (X k (t)), k = 1, 2 are the laws of X fc (t), jfe = 1, 2 under P and W M (- A 
a M ) + J Acrjv/ tt(t) dt is a standard Brownian motion stopped at % under P M . Now, arguing 
as in the proof of Proposition HI we conclude that there exists a uniform constant C < oo 
such that 

/•■Acta/ 

(2.16) \X 1 (-Aa M )-X 2 (-Aa M )\<C / \X\t) -X 2 (t)\ +d 1 (C P (X 1 (t)), C P (X 2 (t))) dt 

Jo 
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holds with probability 1 under P M . However, since for any fixed MeN the measures P M 
and P are equivalent, (1 2 . 1 6 [) must hold under P with probability 1 for all MgN. Taking 
the expectation in (I2.16P under P and taking the limit M — > oo, we obtain 

(2.17) E[\X\-)-X\)\]<C [ E[|X 1 (t)-X 2 (t)|] + rf 1 (£p(X 1 (t)),£ P (X 2 (t)))dt 

J o 

in view of Fatou's Lemma and the Monotone Convergence Theorem. Recalling that 

(2.18) d^CpiX 1 ®)^**®)) < E[\X\t) -X 2 (t)\], 

applying Gronwall's Lemma and using the path continuity of X 1 and X 2 , we end up with 
X 1 ^) = X 2 (-) on [0,T] with probability 1 under P. 

Step 3. Lastly, we need to show the assertion (v) in Proposition [5j By the definition 
of the topology on X, the precompactness of the set {p n : n G N} amounts to showing 
that the sequence p™, n e N is tight with respect to the topology of weak convergence on 
Mi(C([0,T],M + )) and that the set {n(p n ) : n G N} is precompact in y = C([0,T], Mi(E+)). 
We prove these two assertions in steps 3A and 3B, respectively. 

Step 3 A. We need to verify that for every e > there exists a precompact set K e C 
M 1 (C , ([0,T],M + )) such that 

(2.19) supP(p n i K e ) < e. 

ngN 

Moreover, by Prokhorov's Theorem, it is sufficient to demonstrate that for each e > we 
can find a set K € satisfying ( I2.19P such that for every L G N there is a precompact set 
K e ' L C C([0,T],R+) with 

(2.20) sup 7(C([0, T], R + )\K e ' L ) < 1/L. 

To this end, for any function / G C([0, T], ]R + ) and any ( G (0, 1], we let c^(/) be the modulus 
of continuity of / evaluated at (, and write ||.||oo for the supremum-norm on C([0, T],R + ). 
Next, we set 

(2.21) K e ' L = {fe C([0,T],R + ) : c 1/L (f) + \\f\U < M(e, L)}, 

(2.22) 7T = { 7 G M 1 (C([0,T],M + )) : 7(C([0, T],M + )\K e,L ) < 1/L, LeN} 

with a large enough constant M(e, L) < oo to be chosen later. With these choices (I2.19P 
reads as 

(2.23) supPfaLGN: |{1 < z < n : c 1/i (X i (-)) + ||^(-)IU > M(e, L)}\ >n/L) < e, 

nGN V ' 

where X\,Xi, . . . ,X n is a weak solution of the controlled system (12. 5p . Moreover, (I2.23P 
would follow if we can show the stronger bound 

oo 

(2.24) sup VP(c 1/L (X(-)) + ||X(-)||oo > M(e,L)n/L) < e, 

where we wrote £(■) = (Xi(-), X 2 (-), . . . ,X„(-)), and used the notations c\/l and || ■ ||oo for 
the modulus of continuity at 1/L and the supremum-norm, respectively, of a function in 
C([0,T], (M+) n ) with respect to the 1-norm on (M + ) n . 
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We proceed by obtaining an upper bound on E[||^(-)|| 00 ]. To this end, we note that by 
the Cauchy-Schwarz inequality 



X i (-)<X i (0)+ / \b(X l (s),p n (s))\ds+\\g\ 



Xi{s) ds 
+ 



1/2 



Ui(s) 2 ds 

<T(Xi(s))dWi(s) 



1/2 



Taking the supremum over a time interval [0, •] and the expectation on both sides, applying 
the Cauchy-Schwarz inequality again and writing rj(-) := E[sup se r j Xi(s)], we get 



ri(-)<Xi(0)+m sup / \b(X t (0),p n (0))\+C\X l (a)-X l (0)\+Cd 1 (p n (a),p n (0))da 



+ IM 



o 



1/2 



nfslds] kJ 2 + E 



sup 

0<s<- 



o 



with a constant C < oo depending only on the Lipschitz constant of b and 



(2.25) 



k,- = E 



Uiitfdt 



o 



i = 1,2, ... ,n. 



Summing over i, applying the elementary inequality 2a\a<i < a\ + a 2 ,, a±, a<i G R and setting 
r (') = Zir=i r i(') ; we obtain 



r(-) < ^ (X,(0) + |6(X i (0),p"(0))|T+ |MU^/2 + E 

i=l 



sup 

0<s<- 



a{X i {a)) dWi(a) 
+C / r(s)ds 



where C < oo now stands for a constant depending on the Lipschitz constant of b and \\g\ 
only At this point, Gronwall's Lemma yields the estimate 



(t) < ( X M + \KX i (0),p n (0))\T+ WgWooKi/2+E 



i=l 



sup 


I 


-0<s<t 





a(A ? i (a))d^(a) 



£ G [0,T]. Next, applying the L 2 version of Doob's maximal inequality for nonnegative 
submartingales and Ito's isometry, we obtain 



(*) < (^(0) + IK^(0),p n (0))|T+||( ? || oo ^/2+||( 7 || oo 2E 



XAs) ds 



1/2 



e ct , t G [0,T]. 



The elementary inequality 2a < 1 + a 2 , a G R now gives 



r(t) <^(x t (0) + |6(X,(0),p"(0))|T+||^|| oo ^2+||^|| oo )e c *+||^|| oo / r(s)dse ct , te [0,T]. 
i=i ^° 

Applying Gronwall's Lemma once again and recalling (12.91) . we conclude that there exists a 
constant C < oo depending only on b, g, T and the supremum in (12.91) such that 



(2.26) 



n\m-)\U = r(T)<Cn. 
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Next, we give an upper bound on E[ci/l (£(•))]. To this end, we note first that 
J2\Mt) - Xi(s)\ <Y1 / \HX i (a),p n (a))\da + -J2 / [giM^Ma) + u t {af] da 

i=l i=l ^ s i=l ^ s 

„-_i J s 



n , rt 

< s < t < T, 



8=1 



where we have used the elementary inequality 2a\a2 < af + a|, ai,a2 G K. Taking the 
supremum over all < s < t < T with t — s < 1/L and then the expected value on both 
sides and using the fact that b is Lipschitz (see Assumption [1]), we get 

E[c 1/L (X(-))] < i X) IK^(O), P"(0))| + £ r(T) + i X> 

i=l i=l 



2VE sup sup / (j(Xi(a)) dWi{a) 

~l \- k=0,...,[LT]te[k/L,(k+l)/L] 1 J k/L 



with a constant C < oo depending only on the Lipschitz constant of b and ||g||oo- Next, 
applying the L 2 -version of Doob's maximal inequality for nonnegative submartingales, we 
can bound E[ci/l(£(-))] further by 



i=l i=l i=l 



n(T) 



1/2 



< y E 16(^(0)^(0))! + y r(T) + ^f> + ([Lr] + 1)C (n + r(T)), 



«=i t=i v 

where we have increased the value of C if necessary. In view of (12. 9p and (12.261) . this implies 
the existence of a constant C < oo depending only on b, g, T, L and the supremum in (12. 9p 
such that 

(2.27) E[c 1/L (X(-))} < Cn. 

At this point, the bound ( 12.241) with M(e, L) large enough is a direct consequence of Markov's 
inequality, (|2T26|) and d2^7j) . 

Step SB. It remains to check the precompactness of the set {^(p 71 ) : n G N} in y. We 
first show the precompactness of the latter set in C([0, T], Mi(E + )) with the topology of 
uniform convergence, where Mi(M. + ) is endowed with a metric metrizing the topology of 
weak convergence of probability measures. Letting 

(2.28) T — If G C°°(R+) : limsup < oo, ||f |U < oo, \\f" |U < oo}, 

and arguing as in the proof of Lemma 1.3 in [7j, we conclude that precompactness in 
C([0, T], Mi(M + )) would follow if we can show the tightness of the sequence 



(2.29) (P n (-),f):= [ fdp n (-), neN 
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on C([0,T],M) for every fixed / 6 J. In order to prove the latter, we aim to apply the 
criterion in Theorem 1.3.2 of HU and therefore need to show that 

= 0. 



(2.30) V/GJ,A>0: lim lim sup p( sup \(p n (t)J)-(p n (s)J)\>A 

40 n->oo V 0<s<t<T,t-s<e 

Applying Ito's formula, using the union bound and recalling that b is Lipschitz and /', /" 
are bounded for any fixed / G J, we conclude that it is enough to verify 

VAi>0: limlimsupP(-V sup / X^a) da > Ai ) = 0, 

40 n-Hx> \n 0<s<t<T,t-s<eJ s ' 

VA 2 >0: lim lim sup P ( - sup / a(Xi(a)) \ui(a)\ da > A 2 ) = 0, 

40 n^oo V71 0<s<t<T,t-s<e J s ' 

i n r-t 

VA 3 >0: limlimsupPf sup \-T2 /'(A ? i (a))a(A > i (a)) dW*(a) > A 3 ) = 

40 n->oo V 0<s<t<T,t-s<e ' Tl ~^ J s ' 

To prove the first assertion, we apply Markov's inequality in combination with (12.261) to 
bound the prelimit expression from above by Taking limits, one ends up with the first 



assertion. 

To show the second assertion, we first use Markov's inequality and then repeatedly the 
Cauchy-Schwarz inequality to bound the prelimit expression from above by 

— -Vlf sup ( Xi(a)da) ( Ui (a) 2 da) 

WA 2 ^ \-0<s<t<T,t-s<e\ J s / \Js /J 

i=l i=l i=l 

where C < oo is a constant depending only on ||g||oo- In view of (I2.26P and (I2.9p . this readily 
yields the second assertion when one takes limits. 

To deduce the third assertion, we introduce the stopping times 

1 n 

(2.31) tl = inf It > : -^J^(t) > l\, LeN. 

i=i 

Then, the prelimit expression in the third assertion cannot exceed 

p( sup / , (X i (a)) ( r(X i (a))dVr i (a) > A 3 ) + P(r L < T). 

^ 0<s<t<T,t-s<e I U JsAf L J 

We now cover the interval [0, T] by intervals [0, 2e], [e, 3e], [2e, 4e], . . ., use the union bound 
and then apply Bernstein's inequality in the form of Exercise 3.16 on page 153 in [12] (note 
that Problem 3.4.7 in [10] can be applied to the sum of stochastic integrals in the last display) 
to obtain the upper bound 

2T / A\n\ n , 

— exp^-c-^J +P(f L <T), 

where c > is a constant depending only on ||/'||oo and Hfliloo- By taking limits we deduce 
that the left-hand side in the third assertion is bounded above by P(fi < T) for every LeN. 
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Finally, combining Markov's inequality with (I2.26P and passing to the limit L — > oo, we end 
up with the third assertion. 

All in all, we have shown that the set {vr(p n ) : n e N} is precompact in C([0, T], Mi(R+)). 
In addition, we observe that the identity function on M + belongs to J 7 , so that the argument 
above implies that the sequence (p n (-),x), n G N is tight on C([0, T], M). Putting these 
two facts together we conclude that it is possible to extract a subsequence of 7r(p n ), n e N, 
which converges in C([0, T], Mi(R + )) and such that the corresponding functions (p n (-),x) 
converge in C([0, T],R). By Theorem 6.9 in [TO] such a subsequence converges in y = 



In this section we explain briefly how the stochastic differential equations in (II. ip arise 
in the study of stability of the short term interbank lending system. We shall view the 
diffusion (II. ip as a dynamic system that describes the monetary flows between banks from 
the point of view of a financial regulator. Then, taking a particular functional form of the 
drift and diffusion coefficients b, a and letting the number of banks go to infinity, we analyze 
the limiting object of Corollary |2j viewed as the evolution of the distribution of money in a 
large interbank lending system. 

Suppose that Xi(t) > represents the amount of cash and liquid assets (monetary reserve) 
of bank i at time t in a financial system. In order to prevent a financial crisis in the system, a 
financial regulator typically requires banks to hold a sufficient amount of monetary reserves. 
Due to daily banking activities some banks sometimes need additional money to fulfill the 
requirement, while other banks have excess monetary reserves allowing them to lend money 
to others. The interbank lending market exists as an institution, which matches the demand 
and supply of such reserves allowing a bank i to borrow from or lend money to other banks 
j ^ i for a short term (for example, overnight). 

From the perspective of a financial regulator the empirical measure p n {-) = - Y^i=i ^Q(-) * s 
an important quantity to be monitored over time. When the empirical measure is close to the 
point mass 5q, the regulator wants to intervene in the system, whereas otherwise the system 
evolves free of interventions. The day-to-day transactions from one bank to another in the 
interbank lending market are typically quite large, which justifies the usage of macroscopic 
diffusions as in (II. ip as models of the monetary reserves (Xi(-), X2(-), . . . ,X n (-)). In (II. ip 
the drift coefficient b, which stands for the monetary flow to or from bank i, is a function of 
the current monetary reserve Xi(-) and the current empirical measure p n {-), which represents 
the state of the whole interbank lending system; the diffusion coefficient a is a function of 
the current monetary reserve Xi(-) only and stands for the volatility of the monetary reserve 
of bank i due to the daily business of that bank. 

For example, we may consider the following choice of b, a: 



where 8 is a positive constant and <p is a function from Mi(R + ) to M + . The constant 5 
can be interpreted as the (common) growth rate of the monetary reserves and tp(p n (-)) as 
the intensity of the monetary flows between banks given the current state of the interbank 



C([0,T],Mi(R + )). 



□ 



3. A Case Study 



(3.1) 
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lending market p n (-)- The diffusion coefficient a in (13. ip implies that the variance rate of 
Xi(-) is assumed to be proportional to the size of Xj(-). The resulting model reads 

1 n 

(3.2) dXi(t) =(S + -J2 - Mt)) V(p n (t))) dt + 2 ^fxtf) dWi(t), i = 1, 2, . . . , n 



with (Wi(-), W^C"); • • • > W n (-)) being an n-dimensional standard Brownian motion as before. 

Assumption 2. Let tp be a bounded Lipschitz function from Mi(R + ) to M.+. In addition, in 
order to exclude trivial cases, assume that the limiting initial measure p°°(0) = A G Mi(M + ) 
has a strictly positive first moment 



(3.3) 



V dX(y) > 0. 



Note that under Assumption [2] the specification (13. ip satisfies Assumption [T] (with 0(x) 
In this particular case, the sum process 

(3.4) S n (-) :=X 1 (-)+X 2 (-) + ... + X n (-) 

is a squared Bessel process of dimension nS: 



(3.5) 

where 

(3.6) 



S n (t) = S n (0) +ndt + 2 v / ^)d/3(s) 

Jo 



t > 0, 



n „. 



-^(t)\V2 



dWi(t) 



is a standard Brownian motion by Levy's characterization theorem. In particular, this allows 
us to conclude E[S" n (t)] = S n (0) + n8t, t > 0. Therefore, the computation 



E 



^n(T) 



E 



v -X;x,(o)-m A x ' 

i=l 



E 



2 



V^>)d/3( M ; 



+2E 



S n (0) 



n 



m x )\- 

n 



AMd/3(«) 



the convergence of the initial conditions, Ito's isometry and the Cauchy-Schwarz inequality 
show that 



(3.7) 



lim E 

n— >oo 



Sn(T) 

- oT — m\ 



n 



0. 



Hence, by the L 2 version of Doob's maximum inequality for continuous martingales 



Ve>0: P( sup 

■ 0<t<T 



xdp n (t)(x)-St-m x 



>e)< 



-E 



Sn(T) 

oT — mx 



n 



-> 



Thus, the sequence of (random) functions t h- > L x dp n (t)(x), n G N converges in the space 
C([0, T], to 1 1 — y St -\- m\ in probability in the limit n — > oo . Putting this together with 
Theorem [T] and Corollary [2] we obtain the following. 
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Corollary 6. Suppose that Assumption® holds. Then, the X -valued sequence p n , n G N of 
empirical measures corresponding to the particle systems in (13 HJ) satisfies a large deviations 
principle on X with scale n and a good rate function. Moreover, that sequence converges in 
distribution to the law of the unique strong solution of the non-local SDE 



with initial distribution X, where W(-) is a standard Brownian motion and C(X(t)) is the 
law of X(t). Lastly, the solution of ( 13. 81) satisfies K[X(t)} = m\ + 8t, t G [0, T]. 

Note that the SDE (13. 8ft includes as special cases the SDEs satisfied by squared radial 
Ornstein-Uhlenbeck (OU) processes and the SDEs satisfied by squared Bessel processes. 
Indeed, if one lets ip be a constant function and takes 5 = 0, the solution of ( 13. 8p becomes 
a squared radial OU process (see [5], [T7] and [5] for the definition and properties of the 
latter). On the other hand, with the choice p(-) = the solution of (I3.8P becomes a squared 
Bessel process of dimension 8. 

3.1. Non-local squared radial Ornstein-Uhlenbeck processes. In the following, we 
shall consider the effect of interactions through the function ip(-) on the properties of the 
nonlocal SDE (13. 8p when 8 = 0. To wit, 



and C(X(0)) = X. This should be thought of as an approximation of the case when the 
growth rate 8 is small. We will refer to ( 13.91) as the reduced form of (13. 8p . The existence 
and uniqueness of the strong solution to the reduced form SDE (I3.9P can be shown as in the 
previous section (see also [13]). As for the squared radial OU process, the path properties 
of the process X(-) in ( 13. 9p are determined by the bounds of tp. 

Proposition 7. Suppose that Assumption® holds and that X is a point mass at some xq > 
(so that X(0) = xq and m\ = xq). Then, the solution X(-) of (13. 9p has the following 
properties. 

• If m\XD.i fl £M 1 ( p(f J >) > 2, then X(-) is transient; if nix sup^g^ <£>(/i) < 2, then X(-) is 
recurrent. 

• If mx'mi^^ tp(fi) > 2, then X(-) almost surely never hits the origin. 

• If vii\ inf <p{n) > 1 and mASup^^ y?(/u) < 2, then X(-) hits the origin in finite time 
with some positive probability: 



(3.8) 



dX(t) =[8+ (m A + 8t- X(t))ip(C(X(t)))]dt + 2y/X(jjdW{t) 



(3.9) 



dX(t) = (m A - X(t))ip(£(X(t)))dt + 2y/X(t) dW(t), t > 



VT>0: P(X{t) = for some t G [0, T]) > 0. 



However, with probability one the semimartingale local time 




accumulated by X(-) at the origin vanishes identically. 

• If m\vai </?(a0 > and m\ sup^g^ <p(fi) < I, the process X(-) does accumulate local 
time at the origin and the origin is instantaneously reflecting. 
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Proof. We shall transform the process X(-) to a more tractable one by introducing a suitable 
integrating factor. Define 

(3.10) C(0 •-= eotp (| jT *>(£(*(*))) df), ^(-):= j^(t) 2 dt and £(-) == *("K(-) 2 • 

Note that under Assumption [2] the function tp is bounded above by some constant C < oo, 
so that 1 < < exp(Ct/2), t > 0. Hence, ^(-) is strictly increasing, satisfies i < tp(t) < 
J 'exp(Cs) ds, t > and has an inverse In addition, by Ito's formula 

n d£(f) = d(X(t)C(t) 2 ) = X(t) d(c(t) 2 ) + C(t) 2 dX(t) 

= m A p(cm/cm) at) 2 & + 2 yW) cw dww. 

Equivalents, setting f(0) = ^(0)^ = x , D(0^= (^"HOVC^'H-)) 2 )) and intro- 

ducing another Brownian motion W(-) with W / (?/'(-)) = f ((t) dW(t), we obtain 

e(") = ^(0)+ / m A ^(£(^- 1 W)/C(^- 1 W) 2 ))dt+ / 2 V / ^)dW(t) 

= f(0) + / D(t)dt + / 2 v / f(t)dW(t). 

Comparing the latter SDE with the stochastic differential equation for a squared Bessel 
process of dimension d: 

(3.13) dR(t) = ddt + 2^/R(t)dW(t) 1 

we see that the process £ can be intuitively thought of as a squared Bessel process of variable 
dimension D(-) running with respect to the clock ip(-). On a formal level, the statements 
in the proposition can be established by combining the comparison theorem of Ikeda and 
Watanabe (see [9]) and the corresponding properties of squared Bessel processes (see e.g. 
Chapter XI in [12]). For example, if m\ sup <p(fi) <2 , then 

D(-) < m\ sup <p(fi) =. S + < 2 . 

Hence, by virtue of the comparison theorem of Ikeda and Watanabe (see [9]) the process 

£(•):= e(0)+/ S + dt+ / 2^/at)dW(s) 
Jo Jo 

dominates the process £ of (13.121) ; that is, £(■) < £(•) with probability 1. Therefore, since 
for any T > the process £ reaches the origin before time T with positive probability, the 
same is true for £. Moreover, the recurrence of £ follows from the recurrence of £. The other 
statements in the proposition can be shown in a similar fashion. □ 

Below, we shall analyze the dependence of the moments of X(-) on the function tp. There- 
fore, as a preliminary we need to establish the finiteness of such moments. 

Proposition 8. Suppose that Assumption^ holds and that X is a point mass at some x$ > 
(so that X(0) = xq and m\ = xq). Then, the solution X(-) of (13.91) has finite moments of 
all positive orders p > and, moreover, E[(sup 0<s<t X(s)) p ] < oo for all t. Furthermore, 
the stochastic integral J Q X(t) p dj3(t) with respect to any standard Brownian motion /?(•) is a 
martingale of class DL. 
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Proof. We shall use the identity ( 13. lip from the proof of Proposition [7| to estimate the 
quantity E[(sup 0<s<t X(s)) p ] from above for any fixed p > and t > 0. By Assumption [2] 
there exists a constant C > such that sup 6 ^ <p(fJ>) < C- Recall also that we may view 
£(•) intuitively as a time-changed, squared Bessel process of variable dimension 

d(-) = mMmrW/ar'i-)?) < \™ x c] 

running according to the clock ip(-) < J Q exp(Ct) dt =: "0(")- Here, \x] is the smallest integer 
not less than x. Applying the comparison theorem of Ikeda and Watanabe (see [9]) as in the 
proof of Proposition [7J we conclude that the squared Bessel process t given by the strong 
solution of 



(3.14) 

satisfies 
(3.15) 

In particular, 



t(')=£(°)+ / \m x C]dt + / 2y/itfjdW(t) 



sup f (s) 

0<s<t 



e(-)<^(0) and < t(- 



sup £(V> 1 ( , 0(s))) < sup t(s) < sup t(s). 

0<s<< 0<s<V(t) 0<s<^(t) 



Recall now that the running maximum of a reflected Brownian motion on R + has finite 
moments of all positive orders. Moreover, the moments of the running maximum of a 
squared Bessel process of dimension [m^C] can be bounded above by the corresponding 
moments of the sum of squares of running maxima of [m^C] reflected Brownian motions. 
Since X(-) = C(-)~ 2 £(0> it follows that 



E 



(3.16) 



( sup X(s)Y 

L 0<s<t 



= E 
< E 



( sup C(s) 

0<s<t 

( sup £(s)) P 

L 0<s<t 



2 mr 



< E 



( sup t(s)) 1 

0<s<$(t) 



< oo, t > 0. 



Given a fixed a > 0, let us denote by T a the set of stopping times r on the underlying 
probability space, for which P(r < a) = 1. A stochastic integral of the form J (X(t)) p d/3(t) 
is a local martingale and, for every r G 7^, it holds 



E 



X{tf d(3(t) 



< E 



X(t) 2p dt 



< E 



X(t) 2p dt < aE[( sup X{t)) 2p ] 

0<t<a 



< aE 



{ sup t(i)) 

0<t<^(a) 



2p 



< OO 



where the last inequality is due to (I3.16p . Thus, 



sup E 



X(t) p df3(t) 



< oo. 



That is, the family 



[ J\{t) p d(5{t) : reT a ) 



is uniformly integrable. Since this is true for every a > 0, the local martingale j X(t) p d/3(t) 
is of class DL and therefore a martingale (see Problem 1.5.19 in [10]) for every fixed p > 0. □ 
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3.2. Moments. We shall now analyze the dependence of the moments of the solution X(-) 
of (13. 9p on the function ip in the situation of Propositions [7] and [HI An application of the 
last statement of Proposition [HJ yields 



E[X(t + h)- X(t)] = E 



t+h 



(m x -X(s)) V (C(X(s)))ds 



+ E 



t+h 



2y/X(s) dW(s) 



'-Jt 



(m x -E[X(s)})if(C(X(s)))ds 



for all t > 0, h ^ such that t + h > 0. Dividing both sides by h and taking the limit h — > 0, 
one obtains an ODE for the first moment E[X(t)] : 

A E[X(t)} = (m A - E[X(t)]) <p(C(X(t))), t > 

with the initial value E[X(0)] = m x = xq . The solution of this ODE is 
(3.17) E[X(t)] =m x = x , t> 0. 

We conclude that one needs to consider higher moments of X(-) to see the effect of the 
interactions through the function tp. 

Applying Ito's formula to X 2 we obtain 



X(t+h) 2 -X(tf 



t+h 



AX(s)-2 <p(£(X(s))) (X(s) 2 -m x X(s)) 



ds- 



t+h 



4I(s) 3/2 dW(s) 



for all t > 0, h ^ such that t + h > 0. Taking the expectation of both sides and using the 
last statement in Proposition [8] with p = 3/2, one computes 

"t+h 



E[X(t + /i) ] — E[X(ty 



4m A - 2^(£(X(s)))(E[X( S ) 2 ] - ml) 



ds. 



Dividing both sides by h and letting h go to zero, we obtain an ODE for the variance 
V{-) := Var(X(-)) := E[X(-) 2 ] - m\ of X{-): 



(3.18) 



dV(t) 
dt 



4m x -2tp{£{X(t)))V{t), t>0 



with initial value V(0) = 0. Thus, the variance V(-) increases (decreases, respectively) 
in time if ip(C(X(-)))V(-) < 2m x (ip(C(X(-)))V(-) > 2m x , respectively). In other words, 
when the interaction through the function ip(-) is relatively small (large, respectively), the 
variance increases (decreases, respectively). Since V(0) = 0, the variance V(-) increases in a 
neighborhood of t = 0. More precisely, if one defines the constant 

t* := inf{t > : <p(£(X(t)))V(t) = 2m A }, 

the variance increases on the interval [0,£*). In addition, setting </>(•) := (p(C(X (■))) we can 
write the solution of the ODE (I3.18P as 



(3.19) 



V(-) 



4m A exp 2 / <f>(s) ds dt 



■ exp 



(pit) dt 



Am x j exp y - 2 j <p(s) ds^j dt. 



The higher order moments of X(-) can be computed in a similar manner. 
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3.3. A linear first-order PDE. Let us introduce the Laplace transform of the solution 
X(-) of O: 

(3.20) U(t,x) := E[exp(-xX(i))] = / e~ xy £(X(t))(dy), (t,x) G (M+) 2 , 

where £(X(t)) is the law of -AT(t) as before. Note that C(X(-)) depends implicitly on the 
choice of the interaction function <p : Mi(R + ) — > M + . 

In a fashion similar to the arguments in section 13.2} we may apply Ito's formula to 
exp(— xX(t)), integrate over the interval (t,t + h), take the expectation on both sides and 
pass to the limit h J, to obtain 

<9/7 

w (t, X) = ( - M *) + E[X( t )e-<»]) x^ Wi)) ) + 2^ E [X( t )e-<"] , 

(t,x) G (ffi+) 2 . Noting that dU(t,x)/dx = — $L[X(t) exp(— xX(t))] and recalling the notation 
</>(•) = (p(C(X (•))) we can rewrite the latter PDE in form of a linear first-order equation: 

dU dll 

(3.21) —(t,x)+x((/>(t) + 2x)—(t,x)+mxx<t>(t)U(t,x) = 0, (t,x)G(M+) 2 
at ox 

with initial value U(0,x) = f R e~ xy \(dy). 

3.4. Stationary distribution. Suppose that the solution X(-) of (I3.9P is recurrent, so that 
it has a unique stationary distribution a. In particular, this is the case under the condition 
m x sup^;^^ <£(a0 < 2 (see Proposition [7]). In this case, it follows from (I3.2ip that the 
Laplace transform u(x) := J R+ e~ xy a(dy) satisfies 

du 

Up* + 2x) — (x) + ip*m x u(x) = 0, x G R+, 
ax 

where <p* := (f{a). Solving this ODE with the boundary value u(0) = 1 we obtain the 
Laplace transform 



u(x) = {l + (2/^)x)-** mxl \ xe 



of the Gamma distribution with parameters a := 2/cp*, b := <p* m\/2: 
(3.22) a(d!/) = _J_(^ 1 - 1 e Xp (_^) d!/ , y€R+ . 

In addition, the above derivation shows that the function tp : Mi(R + ) — > R + and the 
stationary distribution a of f)3.22p must be coupled through (p(a) = cp*. 
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